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^' ) | Abstract. The OSE (Offline Simulations Environment) simulator of the LPF (LISA 

Pathfinder) mission is intended to simulate the different experiments to be carried out in 
flight. Amongst these, the thermal diagnostics experiments are intended to relate thermal 
disturbances and interferometer readouts, thereby allowing the subtraction of thermally induced 
interferences from the interferometer channels. In this paper we report on the modelling of these 
simulated experiments, including the parametrisation of different thermal effects (radiation 
pressure effect, radiometer effect) that will appear in the Inert ial Sensor environment of the 
LTP (LISA Technology Package). We report as well how these experiments are going to be 
implemented in the LTPDA toolbox, which is a dedicated tool for LPF data analysis that will 
allow full traceability and reproducibility of the analysis thanks to complete recording of the 
processes. 



1. Introduction 

The Thermal Diagnostic Subsystem of the LISA Technology Package (LTP) onboard LISA 
Pathfinder (LPF) [H [2] includes the generation of thermal disturbances in some thermal sen- 
sitive locations of the LTP to induce perturbing signals in the interferometer (IFO) readouts. 
The objective of such controlled perturbations is to obtain information of the effects that po- 
tential temperature gradients in the LTP may induce in the final IFO readout during the main 
operational phase, enabling the possibility of removing thermally induced noise from the signal. 

The LISA Pathfinder main noise requirements in the LTP state for the test masses residual 
differential acceleration: 



^lpfH^xIO" 14 



1 + 



lo/2tt \ 2 



3mHz 



ms- 2 Hz- 1/2 (1) 



in the bandwidth between 1 mHz and 30 mHz [3 J 



Therefore, a set of Thermal Diagnostics experiments consisting of the injection of controlled 
loads of heat at specific spots must be carefully designed, producing in-band perturbing signals 
with high Signal-to-Noise Ratio (SNR « 50). Moreover, a model linking both the heat 
propagation and temperature distributions to final IFO readouts must be designed, and it is 
intended to be introduced into two LPF key simulators: 

• The OSE (Offline Simulation Environment) simulator: Located at ESAC (European Space 
Astronomy Centre), this simulator is intended to validate all LPF mission operations 
considering the different interfaces and data transfer processes. 

• The LTPDA simulator: Developed by the LTP Data Analysis (LTPDA) team, within the 
LTPDA software [3j. It is a linear simulator based on State-Space Modelling (SSM) intended 
to analyse the data delivered by the mission. 

Part of this model is presented here in this paper in Section [3j and some first results are 
reported in Section [U Nevertheless, it is important to first introduce the different Thermal 
Diagnostics experiments and the architecture of the Data Analysis associated to them. 

1.1. Thermal Diagnostics Experiments 

The LTP is equipped with a set of 14 heaters and 24 temperature sensors (see Figures Q] and [2j) . 
Different experiments are intended to apply heat signals through the heaters and record 
the temperature readouts. These experiments are gathered in the Experiment Master Plan 
(EMP) [3] which aims to characterise the thermal effects in Table Q] by switching on different 
heaters at each case, and consequently thermally disturbing different parts of the LTP. 



Table 1. Types of Thermal Diagnostics Experiments. Heaters, associated thermal effects and 
expected main consequences. 



Heaters activated 


Thermal effects 


Main consequence 


Electrode Housing (EH) 


Radiometer effect, 
radiation pressure 


Force and torques 
on Test Masses (TMs) 


Optical Window (OW) 


Glass and clamp 
thermal expansion 


Optical Window glass 
refractive index change 


Suspension Struts (SS) 


Optical Bench asymmetric 
thermal expansion 


Optical Bench (OB) 
thermoelastic stress 



1.2. Thermal Diagnostics Data Analysis 

A complete model simulating the different consequences of the experiments (temperature 
responses, dynamics on the TM, and phase shift in the IFO) on one hand, and the telemetry data 
of the IFO and the temperature sensor readouts, on the other, will be used in the subsequent 
data analysis, when all the data will be collected and processed so more accurate values can be 
obtained for the model coefficients. 

The global scheme is presented in Figure [3l where blue blocks refer to the thermal 
model-obtained temperature, while red ones refer to the blocks representing conversion from 
temperature distributions to the different perturbations that can disturb the readouts of the 
interferometer. 
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Figure 1. Scheme of the Thermal DS (Diagnostics Subsystem) items locations in the Electrode 
Housings, on the Optical Windows and on the Optical Bench, with their enumeration. 




Figure 2. Scheme of the Thermal DS items location on the Suspension Struts. CAD rendering 
of the LTP courtesy of EADS Astrium GmbH. 



The Data Analysis work must handle two different types of simulations: (1) a simulation 
providing data that will be observed during the mission, as temperature sensor readouts, and 
(2) a simulation providing temperatures at non-observable points whose temperatures are inputs 
to the blocks converting temperature fluctuations into IFO perturbations. Block (1) will be used 
to validate the thermal model used in (2). 

The modelling of the different blocks in Figure [3] varies depending on the feasibility of carrying 
out experimental on-ground campaigns able to reproduce the disturbing thermal effects of each 
experiment in Table [TJ Consequently, experiments for thermal qualification purposes on the 
Optical Window and on the Optical Bench were planned, the first carried out in AEI Hannover 
Laboratory (2007, [6J) and the latter taking place at the time of writing in the frame of the 
LPF Thermo-Optical Qualification Model (TOQM) campaign in IABG facilities in Ottobrunn, 
near Munich, being run by EADS Astrium- Germany and EADS Astrium-UK (November 2011, 
[7J). As thermal effects on the free falling TMs can not be directly tested, an accurate model for 
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Figure 3. Scheme of the Thermal Diagnostics model. LPF picture courtesy of ESA. 



the EH thermal effects is being developed, the results of which are going to be compared with 
similar experiments carried out at Trento University with a torsion pendulum [8j. 

2. Thermal effects in the EH 

Temperature differences across the TMs environment can produce differential pressures that turn 
into net forces and torques on the TMs. Three different thermal effects have been identified [H] 
as possible mechanisms that will create noticeable dynamic effects on the TMs: 

2.1. Radiation pressure effect 

The radiation pressure effect is based on the temperature-dependence of the radiation emitted 
by a surface, following the Stefan-Boltzmann Law. The force applied on a surface j by the 
presence of another surface i at a different temperature is expressed as: 
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where a is the Stefan-Boltzmann constant, c is the speed of light, S pp is a system geometric 
factor from surface % to surface j which contains information of both the system geometry and 
potential multi-reflections, q is the emissivity of surface z, Aj is the area of the section where the 
force is applied, Tq is the absolute temperature and ATij is the temperature difference between 
two surfaces, written as ATij — T{ — Tj. 

2.2. Radiometer effect 

The radiometer effect appears in rarefied atmospheres where the particles have a mean free path 
much longer than the distance between the surfaces of the enclosing volume, i.e in systems with 



large Knudsen number. The consequent force is represented here as: 

F^ = \^il ATzj [N] (3 ) 

where p stands for the remaining gas static pressure in the enclosure and S^ M is a system 
geometry coefficient for the radiometer effect, which includes multi-reflections as well. Aj, Tq 
and ATij have the same meanings as in Equation [21 

2.3. Outgassing effect 

The outgassing effect is caused by the detachment of gas molecules from the surface, in a 
low pressure environment. Its behaviour is strongly related with the surface roughness and 
the material properties; although its consequences could be similar to the ones from the other 
thermal effects, it is not yet included in the model because of the uncertainty of some of its 
parameters. 

More details of all these effects can be found in [H] and the values considered for their main 
parameters are presented in Table [21 where the radiative emissivities presented there refer to 
the TM and EH surfaces Q3J]. 



Table 2. Default values considered for the thermal effects expressions. 



Symbol 


Parameter 


Value 


P 


Environment static pressure 


10~ 5 Pa 


To 


System temperature 


293 K 




TM surface radiative emissivity 


0.02 




EH surface radiative emissivity 


0.05 



3. Thermal effects modelling 

In this Section, we focus on how the thermal effects inside the Electrode Housing are being 
modelled. Such a model can be divided into two main blocks, a thermal block that provides 
temperatures at specific spots when heaters are activated and another block that implements 
the different thermal effect expressions presented in Section [21 taking into account the geometry 
of the system and other effects as the potential multi-reflections between surfaces. 

3.1. ESATAN model block 

A thermal model of the LTP in ESATAN software (developed by Carlo Gavazzi Space) together 
with a general model of the whole spacecraft is maintained at ESTEC (European Space Research 
and Technology Centre). It provides node-to-node transfer function samples, where the input 
is the heat power applied to a node and the output is the temperature response at another 
requested node. Temperature responses are calculated as 

T s (s) = H h ^ s (s) ■ q h (s) (4) 

where is the heat power applied to heater /i, the transfer function and T s the 

temperature response at sensor s, in Laplace Domain. 

Different representative nodes in the thermal model were selected, and series of transfer 
functions relating all the DS heaters to these interesting nodes were obtained and implemented 
into a single state-space model. 



3.2. Geometry and multi-reflections considerations 

Since the real system is not composed of infinite parallel surfaces, it is necessary to consider 
finite size effects when modelling the EH. 

First of all, the surfaces of the system composed of the TM and the inner walls of the EH need 
to be discretised properly so potential temperature gradients across a surface can be represented. 
The discretisation adopted is shown in Figure [H where each TM/EH surface is divided into 4 
new facets, thereby enabling the representation of both forces and torques in the three space 
axis. 
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Figure 4. Discretisation considered for the temperature distribution in the inner faces of the 
EH (exploded view) and on the TM with numeration of the facets. Each facet is considered 
isothermal. 



On the other hand, in order to model radiation reflection effects in the EH it is necessary to 
develop a multi-reflection model for our specific system. Additionally, the model must include 
the possibility of generating not only pressures but also shear forces on the TM surfaces, as they 
are also important [8j. 

Such conditions have required the calculation of all facet-to-facet force contributions coming 
from the different surfaces in a vectorial way (through the computation of the Radiative 
Exchange Factors [H]), so the different resulting components can be disentangled on each 
surface. The main trade-off of these considerations is the assumption that all reflections are 
diffuse [12]. 

This multi-reflection model is then extrapolated and applied to the Radiometer effect, where 
equivalent Geometric Exchange Factors have been considered in replacement of the Radiative 
Exchange Factors, which only apply to the Radiation Pressure Effect. The parameters used in 
such extrapolations are still under study, but first results are quite satisfactory (see Section Sj). 
Table El presents the current values considered for the different effects and surfaces. 

Finally, as its implementation is linear and it has no dynamics associated, this block can be 
considered a post-processing part of the temperature state-space model block. 

4. Results 

Under the assumption of constant temperature at the TM, we report here first results obtained 
from the implementation of the model. Specifically, we present forces and torques obtained from 
the simulation of one planned run from the EMP [13j. In this case, the exercise is applied to 
EH1 and the heaters activated are HI and H2. The parameters of the input signal are given in 
Table [4| and Figure [5] shows graphically the input signal shape. 



Table 3. Parameters considered in the multi-reflection model. Radiative emissivities have 
already been presented in Table 



Symbol 


Parameter 


Value 


Ptm, rp 


TM radiative reflectivity 


0.98 


Peh, rp 


EH radiative reflectivity 


0.95 


£tm, rm 


TM equivalent emissivity for radiometer effect 


0.02 


£eh, rm 


EH equivalent emissivity for radiometer effect 


0.05 


Ptm, rm 


TM equivalent reflectivity for radiometer effect 


0.8 


Peh, rm 


EH equivalent reflectivity for radiometer effect 


0.8 



Table 4. Input signal parameters. 



Heaters activated 


HI and H2 alternatively 


Period 


2000 s 


Duty cycle 


50.00% 


Power 


10 mW (5mW to each physical heater) 


Duration 


4000 s 


Total signal duration 


20000 s 



0.01 
0.005 
0, 
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Figure 5. Input signals. HI and H2 are activated alternatively. 



Figure [6] shows the evolution of the EH temperature and the temperature differences, as 
measured by the sensors in the indicated faces. Temperature sensor readouts in the same face 
of the EH are not identical, but their differences are negligible in this context. 

Note the following temperature evolution features: 

(i) When the heater HI on the +X face is activated, the temperature of the latter starts to 
rise until the heat input is discontinued, during the second half of the heating signal period. 
At the same time, the -X face also feels the heating, and its temperature increases, though 
significantly less. 
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Figure 6. Plots of temperature increments and differences on EH1 facets. The temperature 
differences have been obtained by previously averaging the temperature at all the facets from 
each wall. The facet position follows the notation in Figured! 



(ii) During the second half-period, when HI is off, the +X face cools down, but does not end 
up at the original temperature. The activation of H2 in the -X face actually causes the 
temperature of the +X face to slightly start rising towards the end of the first H2 cycle. 

(hi) As the process continues, there is a consistent drift upwards in the temperature of both 
faces. When the heaters are completely switched off, the temperatures of the faces match 
each other -approximately from 5000 seconds-, then they continue in parallel a smooth 
decline which asymptotically approaches the initial temperature. 

In Figure El we see as well plots of the temperature differences between the different faces. 
As can be observed, a very symmetric pattern shows up in the case of +X and -X faces. 

Finally, the forces and torques obtained on TM1 considering the contribution of each different 



effect are presented in Figure At each plot the contributions of each thermal effect are 
presented separately, and the parameters applied are the default ones used in the algorithm. 
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Figure 7. Plots of forces and torques on TM1, with the disentangled contributions of the 
different thermal effects under the default parameters of the algorithm. 



5. Conclusions 

A model for the thermal Diagnostics Experiments in the Electrode Housing to be carried out in 
the LPF has been presented. Improved from older versions, this model can now represent multi- 
reflections for the different thermal effects in a vectorial way, and it includes representations of 
the radiation pressure and the radiometer effect. 

First results show that noticeable thermally-created forces can be obtained in the IFO- 
sensitive axis -X axis- and that the rest of perturbations, i.e. forces on the Y and Z axis 



and the different torques, present much smaller effects. On the other hand, all torques observed 
are negligible as they are far below the DC 1.1 x 10 -11 Nm specified in [Hj. Additionally, an 
interesting temperature gradient is observed on the Z axis, though it is out of the bandwidth of 
interest. 

However, some of the parameters considered for the model are still under on- going study so 
improved values are expected to be obtained in the future. These parameters refer mostly to 
the multi-reflections block. 

Future work also includes the integration of all the Thermal Diagnostics Experiments into a 
single Module able to simulate the thermal disturbances in the LTP caused by the different DS 
heaters, i.e. adding to the current model the OW phase shift and potentially thermoelastic effects 
on the OB caused by the OW and the Suspension Strut heaters respectively. The implementation 
of the part involving the EH and the OW in the LTPDA simulator so as to obtain IFO readouts 
from the thermal experiments is already under way. 
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